Exact Solution of the strong coupling t-V model 
with twisted boundary conditions 
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We present the solution of the one-dimensional t-V model with twisted boundary conditions in the 
strong coupling limit, t ^ V and show that this model can be mapped onto the strong coupling 
Hubbard chain threaded by a fictitious flux proportional to the total momentum of the charge 
carriers. The high energy eigenstates are characterized by a factorization of degrees of freedom 
associated with configurations of soliton and antisoliton domains and degrees of freedom associated 
with the movement of "holes" through these domains. The coexistence of solitons and antisolitons 
leads to a strange fiux dependence of the eigenvalues. We illustrate the use of this solution, deriving 
the full frequency dependence of the optical conductivity at half-filling and zero temperature. 



PACS numbers: 

I. INTRODUCTION 

The extended Hubbard model and its spinless version, 
the t-V model, have been extensively studied due to 
their relevance in the comprehension of the behaV|ipr of 
strongly correlated compounds such as cupratestJEl and 
organic conduct or sJj'lI Much of the present understand- 
ing of these models has been a consequence of the ex- 
act solutioffl|-jn one dimension by the Bethe ansatz (BA) 
technique B'LI The evaluation of the correlations remains 
however a hard task within the Bethe ansatz framework. 
For the Hubbard model, further progress was possible 
in the strong, coiupling limit due to the simpler form of 
the solutionJjTiS The eigenfunctions in this limit factor- 
ize as a product of a wavefunction of non-interacting 
spinless-Armions and a wavefunction of a squeezed spin 
chain.EnlS This spin-charge factorization simplifies the 
calculation of correlations and in particular, it has bee: 
used to determine thajxipmentum distribution function, 



the spectral function 
lower Hubbard band, 



the sum rules of the uppet-and 
iil and the Green's functionO of 
this model. An alternative solution tp-that of the Bethe 
ansatz was also possible in this limit .E3i]j 

The t-V model is apparently simpler than the Hub- 
bard model due to the absence of spin degrees of freedom. 
This model can be mapped onto the anisotropic Heisen- 
berg model (more precisely, the XXZ or Heisenberg-Ising 
model) by the Jordan- Wigner transformationjla whose 
Bethe ansatz solution has long been known.EZl In the 
strong coupling limit, the t-V model, despite its appar- 
ent simplicity, remains somewhat foggier than the lEiirb-i 
bard model. For instance, the Bethe ansatz solutionll3"E2l 
presents us eigenvalues expressions with phase terms 
whose physical meaning is not clear. Another curious 
fact is that the Luttinger liquid exponenttJ that char- 
acterizes the low energy excitations-of the strong cou- 
pling t-V model is density dependentE3 in contrast to the. 
strong coupling Hubbard model where it is a constant .□ 
Since this exponent is closely related to Fermi surface 
phase shifts (a holon Fermi surface in the case of the 



Hubbard model) jl it is worthwhile to investigate how 
these phase shifts will be modified. In this paper, we 
present a non-Bethe-ansatz solution for the strong cou- 
pling one-dimensional t-V model which is closely related, 
to the solution of the strong coupling Hubbard modeltil 
and which clarifies the previous issues. The simple factor- 
ized form of this solution (and the low degeneracy of the 
eigenvalues) will, we believe, allow an easy calculation of 
correlations. 

The t-V Hamiltonian for a ring with L sites is 



H 



-t 



(1) 



is the fermion creation operator on site i, 
and V is the nearest-neighbor Coulomb inter- 



where c| 

Til — C^ C{ 

action. The one-dimensional strong coupling t-V model 
(as the Hubbard model) is a classic example of a system 
which exhibits a metal-insulator transition upon doping. 
lit — 0, the fermions are localized and all states with the 
same number of pairs of nearest-neighbor occupied sites, 
^-n^ni+i, are degenerate. This degeneracy is lifted if 
t/V is finite and up to first order in the eigenvalues 
are obtained diagonalizing the Hamiltonian within each 
of the degenerate subspaces. In the strong coupling limit 
i <C T^, we obtain therefore the projected Hamiltonian 

H = - "-i+2)cj+iCi(l - ni_i) + he] 

i 

-i^[ni+24-i-iCi'^i-i + /ic] 

i 

-KF^n,n,+i. (2) 

i 

This limit corresponds to the Jx = Jy Jz limit of the 
anisotropic Heisenberg model. The set of eigenstates and 
eigenvalues of this model can be obtained without having 
to resort to the Bethe ansatz, as we will show below. 

The behavior of the t-V model in the strong coupling 
limit has provided |Siijjoort for a recent conjecture by Zo- 
tos and Prelovsek.E^tJ According to these authors, the 
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t-V model and the Hubbard model at half-filling are per- 
fect insulators, this meaning that the Drudc weight {Dc) 
in the thermodynamic limit remains zero even at finite 
temperature. In particular, in the strong coupling limit, 
they argue that Dc = even for finite systems. Such be- 
havior was confirmed by Peres et al,t3 applying the Bethe 
ansatz method to solve the t-V model in the strong cou- 
pling limit. Here, we present a different solution which 
allows an easy confirmation of the previous results and 
makes clearer the physical picture in this limit. 

The remaining part of this paper is organized in the 
following way. In section II, the low energy eigenstates 
of the model are found. First, we study the one-particle 
problem and present a new path for its solution. Then, 
we show that this solution can be extended to the case 
with N particles. We also find the flux dependence of 
the eigenvalues. In section HI, the general solution is 
presented both for periodic and twisted boundary condi- 
tions (finite flux). We also comment on the higher order 
corrections. In section IV, we compare our results with 
those obtained with the Bethe ansatz technique. The 
transport properties of the model are studied in section 
V. Finally, we conclude in section VI. 

II. LOW ENERGY SUBSPACE 

Let us consider two consecutive sites and therefore, 
nearest neighbors of each other. There are four different 
configurations for this pair of sites, which we will call 
links and they are 

(••);(-•);(•-);(--) (3) 

where a dot stands for an occupied site and a dash for an 
empty one. The total number of these links in the chain 
is equal to the number of sites L, 

N,,+N^, + N,^+N^^=L (4) 

and iV_ , = TV, _ . Further conditions result from count- 
ing the number of holes or particles, 

N,,+N^,=Ne. (5) 

In the limit V/t — > oo, the number of links (••) is 
a conserved quantity and consequently also are , 

and N So, in the strong coupling limit, the model 

merely exchanges the positions of these links. Note that 

A. One particle 

Let us start with the simple case of a single particle 
in our periodic chain. In this case, the interaction term 
is zero and we have a one-particle tight-binding model. 



whose solution is trivial. We are going to solve this model 
in a different fashion, considering as our mobile particle 
the link (— •). This link moves exchanging its position 

with a link ( ). Note that we have one link (— •) and 

L—2 links ( ), and therefore, the total number of these 

links is i = L — 1. 

First, let us define our states in terms of the position 
of this link, 

N>=4J0) (6) 

with 1 < i < L — 1. Note that there is a state c||0) which 
is not included in the previous set of states, but this state 
can be written as 

r-i|i) = T-i4|o) (7) 

where T is the one-site translation operator. The single 
particle Hamiltonian rewritten using this new notation 
becomes 

H/i-t) = E N + 1) (^1 + T-'\i) m + h.c. (8) 

with L = L ~ 1. We now introduce an over-complete set 
of states, constructing from a state a state invariant 
by translation with momentum fc. 

These states diagonalize the Hamiltonian, but we are go- 
ing to proceed as if they did not and as if they were an 
orthogonal set of states. The Hamiltonian becomes 

+e'^\l,k){L,k\ + h.c^ (10) 

The Hamiltonian in a given k subspace has become that 
of a tight-binding model in a chain of i — 1 sites with a 
fictitious flux e^^ . The eigenvalues are given by 

E{k, k) = ~2t cos - (11) 

with k — fi ■ 2t: / L, n = 0, ... L — 1 and k = n ■ 27r/L, 
n — 0, . . . L ^ 1. But clearly, not all combinations of k 
and k correspond to real eigenstates of the Hamiltonian. 
The eigenstates will be of the form 

|fc,A:) = ^Ee''''^^-'|j''^)' (12) 

but obviously, this is a combination of states which are 
the same state except for a phase, that is. 



2 



\j,k} = e''\j + l,k) 



(13) 



and the linear combination will be zero unless the follow- 
ing condition is satisfied 



A;- 



- fc = (mod 27r) 



(14) 



which implies k = kL/L (mod 27r). This condition is 
equivalent to stating that the final state must obviously 
have a momentum k. Note that the previous equation 
has precisely L solutions and therefore, the usual set of 
tight-binding eigenvalues is recovered. 

Let us make a few remarks concerning the above pro- 
cedure. Let {\i)}, i = 1,. . .L, be an orthogonal set of 
states which constitutes a basis for the states of a given 
system and let Hij be the matrix elements of the system 
Hamiltonian between the states \i) and The set of 
cigenstates {lei,;}}, ?' = 1, . . . L, of the Hamiltonian can be 
written in this basis as = Assume now that 

two or more states of this basis were in fact the same 
state. This would lead to a reduced matrix Hij which 
would have the same elements, except for the transitions 
to and from the state which remained from the set of 
equivalent states. These transitions are multiplied by 
the total number of equivalent states. The eigenstates of 
H remain eigenstates of H, but now they constitute an 
over-complete basis of the Hilbert space, that is, they are 
not all linearly independent. They may however remain 
an orthogonal set if the states in excess are identically 
zero as in the case studied above. 



that starting from a state as above, these states appear 
when a link (— •) is at site L, or equivalently, a particle 
is at site L and hops to site 1. However, this hopping 
term can written as 



(17) 



where f is the one-site translation operator in the chain 
of L sites. Note that this translation operator also insures 
that the other pseudo-particles remain in the same sites 
in the reduced chain. Given a state |oi, . . . , ojVe), we 
build as previously the state invariant by translation with 
momentum P, 

1 ^ 

|ai, . . . , aN^,P) = Y.e'''^fi-^\au • • • , a^vj- (18) 

i=i 

The mapped Hamiltonian in the subspace of states with 
momentum P becomes 



^(^) = (^E + ^'''A^ij + (19) 

So, we have mapped the Hamiltonian onto a tight- 
binding chain with L — N sites threaded by a fiux P 
with A^e particles. The eigenvalues are given by 



£;({fc},P) = -2i^cos fc 



i=l 



P 



(20) 



B. several particles 

Let us consider now the case of iVg particles in a chain 
of L sites, but distributed so that there are no links (• •). 
These states are of the form 



|ai,...,ajve) = n(l 



i)4jO) 



(15) 



with {aj} an ordered set of non-consecutive numbers 
(ai > 1). The total number of ( ) and (— •) links 

is L = L — Ne- This state can be mapped onto the fol- 
lowing state of a chain with L — sites, 



|ai,...,a;ve) = n4l0) 



(16) 



i=l 



such that if the first site of this pseudo-chain is empty, 

the first link of the L sites chain is ( ), if it is occupied, 

the first link is (— •). The same reasoning applies to the 
other sites. Note that as in the previous case, there are a 
few states which are not included in the above set, namely 
states where a link (— •) is divided between sites L and 1. 
These states have a particle at site 1 and one should note 



with k = n-2TT/L, and P = n-27r/L, with n = 0, . . . , L — 1 
and n = 0,...,L — 1. Again, not all combinations of 
pseudo- momenta {fc} and P are possible and applying 
the same procedure as before, we arrive to following con- 
dition 



^(fc, --^)-P = (mod27r) 
i=i ^ 



(21) 



which implies 



P^ 
L 



^fc (mod 27r). 



(22) 



The factor L/L converts the total momentum of our 
chain of L sites in the total momentum of the pseudo- 
chain. Note that the set of pseudo- momenta {fc} is not 
enough to define the total momentum P since there may 
be two values of P such that PL/L (mod 2n) is the same. 
In fact, if P = n • 2n/L , with n = 0, . . . , L — 1, the states 
I {fc} , P) with n in the interval [0,L — 1 — L] have the same 
energy as the states |{fc}, P+ ^j^)- In the particular case 
oi L = L/2, given a state with momentum P, one always 
has a state with momentum P + with the same set 
of wavenumbers. The L = L/2 case corresponds to the 
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half filling and indeed one knows that the ground state 
is doubly degenerate, one state having zero momentum, 
the other tt. This is also valid for excited states (with 
N,, ^0), as we shall see in the next section. 

An external magnetic flux (jj can be introduced in the 
problem with the transformation t te**^/^ . The Hamil- 
tonian remains invariant by translation and all the previ- 
ous steps can be repeated, leading to the following mod- 
ification in the eigenvalue expression 



Ei{~k})^E(\~k 



(23) 



The ground state energy is given by Eq. If is 
odd, all single-particle states with pseudo-momentum k 
between ±2-k/L ■ {N^ — l)/2 are occupied and ^ fc = 0. 
Therefore, 



ji^odd 



iVe — -I- 

^ f2TT\ 

-2t > cos I —I 1 

jv.-i \ L J 



sm 



-2t- 



( 



(24) 



If Ne is even, all states with k between — 27r/L- (iVg — 2)/2 
and 2Tr/L-Ne/2 or between -2n/L-Ne/2 and 2Tr/L-{Ne- 
2)/2 are occupied and ^ ^ ztn/L ■ N^/ L. So, 



-2t 



> cos I —i 



L L 



sm 



-2t- 



(^) 



sm 



(i) 



(l 



(25) 



This slight energy difference between thcr-two cases had 
already been pointed out by Kusmartsev.L3 In the pres- 
ence of a small flux 0, the odd expression should be 
multiplied by a factor cos{(j)/L), while for even, a 
—(fi/L term should be summed to the argument of the 
cosine. 

This phase shift between the wavenumbers of the 
ground states with iV and N + 1 particles should be re- 
sponsible for the orthogonality catastrophe in the ther- 
modynamic limit which, for example, leads to a zero 
renormalization constant Z characteristic of a Luttinger 
liquidEiJ (see Ref.0 for a detailed calculation in the case 
of the strong coupling Hubbard model). The renormal- 
ization constant Z is given by the overlap between the 
ground state with -|- 1 particles and the ground state 
with N particles plus a particle at Fermi momentum, 

Z = \{^Gs{N + 1; P = %)|ct JV'G5(iV; P = 0))p (26) 

yielding zero in the thermodynamic limit. Our results 
above indicate that the phase shift depends on the den- 



coupling Hubbard model, the phase shift of the holons 
wavenumbers is independent of the band filling, S — ^. 
That phase shift results from a {—kp) momentum contri- 
bution from the spin sector J3 Here, the phase shift is due 
to the total momentum of the charge carriers. This de- 
pendence on the band filling is in agreement with the fact 
that the anomalojus exponent of this model is indecjiband 
filling dependents The Luttinger liquid velocitiescil that 
characterize the low lying excitations have been found 
by Gomez- SantoaHj for the strong coupling t-V model 
in the thermodynamic limit based on very simple argu- 
ments (basically, the reduction of the effective size of the 
chain) . These velocities and its finite size corrections are 
easily obtained from the previous equations. In the large 
L limit, the Gomez-Santos results are reproduced: 



2t 



VN = — 



L d^Er _ 
n~dW ^ (1 - pf 



vj 



Ld^Er 

TT 9(^)2 



2t{l- p) sin 



2t 



pi: 



piT 



(27) 



(28) 



(29) 



sity, 6 



np. Recall that in the case of the strong 



where P — ^ and vn, vj, and vs are respectilvely the 
particle, current, and sound-wave velocities. 



III. GENERAL SOLUTION 

Let us consider the general case where one may have 
both links (• •) and ( ). First, note that a phase sep- 
arated state (one domain of holes and one domain of 
particles) has no mobile entities in the strong coupling 
limit since any hopping of a single particle would imply 
the breakup of a link (••). So, phase separated states 
will be eigenstates of the strong coupling Hamiltonian 
with eigenvalues given by E = V ■ N,,. Furthermore, the 
same applies to states with several domains if the only 
links (— •) present are the domain walls. Clearly, a hole 
(particle), in order to be able to move, must be within a 
particle (hole) domain. If for a chain with L sites and 

particles, we fix A^, , and A^ , it is the configuration of 

these links that will define how many mobile links (— •) 
one has and consequently, the number of sites L of the 
effective chain for these mobile links. These mobile links 
will move exchanging their position with links (• •) and 
(--)• 

It will prove itself useful to do the following mapping: 

(••) = !!>, 
(--) = U>, 
(-•)-|o), 

with the exception of the links which are domain walls. 
That is, we will map the states of the spinless chain with 
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L sites and Ne particles onto states of a spinful chain 

with L sites and N, , particles with spin up and N 

particles with spin down. The first two links are called 
respectively a soliton and an antisoliton. A general state 
is written as 

|ai,...,aArJ=n4jO) (30) 

1=1 

with {ai} an ordered set of integers between 1 and L. 
Note that now a particle may occupy the first site and a 
link may be divided between sites 1 and L. These states 
will now be mapped onto the states of a reduced chain 
with the number of sites being 

L = iV..+iV__+Ar_.-Arx| (31) 

where N^^ is the total number of domain walls in the 
sequence of spins obtained by the mapping above. The 
above relation leads to the following relation between the 
real and effective chain sizes 

Z = L-7V_.-iVxT (32) 

which reflects the fact that our moving "particles" are 
now the links (— •) with the exception of the ones which 
are domain walls. The two sites that compose such a 
link are effectively reduced to one (or zero, if the link 
is a domain wall), with the consequent reduction of the 
chain effective size. Note that L is always larger than Nf. 
or Nh. 

The state given in Eq. ^ corresponds to the following 
state of the reduced chain 

\ai,...,aN^+N^;cfi,...,aN^+N^)= J| 4,ctJ0> (33) 

1=1 

If site 1 is empty, the first link of the chain of L sites is 
(— •) and in order to have a well defined mapping, we 
impose the condition that first two sites of the chain of L 
sites correspond to the link and therefore the first site is 
empty while the second is occupied. The same applies in 
case of the site 1 being occupied. Links which are domain 
walls are not mapped to the reduced chain (see Fig. 1 
for an example of the mapping). This condition agrees 
with the definition of states of the previous section and 
furthermore, it also implies that certain states are not 
included in the mapping, but, as previously, they can be 
written as translations of states included in the mapping. 
These states which need to be translated appear due to 
hoppings between sites 1 and L, but also sites 1 and 2. As 
previously, we construct states invariant by translation 
with total momentum P, 

1 ^ 

|{a},M,P) = |{«},P) = Y.e^P=f=-'\{a}). (34) 



and keep the same mapping. The states which need to 
be translated lead to e^*^ terms in the mapped Hamil- 
tonian. So that one does not need to be concerned with 
the reordering of operators in the real chain, we will con- 
sider Ne odd. The even case can be solved with minor 
modifications of the procedure below. Let Ng = N^+Np 

Note that in general, the hopping of an electron implies 
simply that aj — > aj ± 1 for some j. Hoppings of a 
particle from 1 to 2 or 1 to L are however more complex 
processes in the reduced chain. In the following tables, 
we describe the action of these hopping terms. In the 
first column of each table, one has the initial state and 
in the last column, the final state after the application of 
the hopping operator. An extra intermediate column is 
present if the final state can not be directly mapped onto 
a state of the reduced chain. The second line in each row 
shows the states in the original chain while the first line 
shows the equivalent states in the reduced chain. 

i) Let us first consider the jump of a link from 1 to L. 
Note that this implies a cj^ci hopping for a link (• •), but 
a. c\cl hopping for a link ( ). 







• • • • • — • — 


— • • • • • • 
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no mapping 






o • • 
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• — 




-1 _ 


• — 
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i • 
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no mapping 






o • • 
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■ ■ • 


• — • • • 


= T 


-1 _ 


- • — 




i- 


• • o o 


no mapping 




e^^ 


o • • 


oi 




• — • — • 


• — 


= T- 


-i _ 


• — 


• • ■ — • — 



ii) Now, the jump of a link from 2 to 1: 



°i ••• i 




— • ... — 


• — — 






— • — • • • • • 


— • — 


i • • • 


i o • • • o 


— • ■ • • — • 







no mapping 


e'^ T ° ■ • • T 










no mapping 


e-'^ T ° • • • i 






= r+U . - . 



The last two cases also occur if the last pseudo-spin is 
not at site L. The Hamiltonian Hi = H — V ■ N,,, in 
the mapped Hilbert space (in the subspace of momentum 
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P), becomes 



where 



_i\Nh l/2(l+ai-aN,)iP 



(40) 



(-1) 



(35) 



with 



hi 
hi 

Hi 



and tio- = t, in the other cases. This is the U = oo Hub- 
bard chain pierced by a magnetic flux. The Hamiltonian 
does not change the sequence of spins {cr}, but circularly 
permutes them. Note that (-1)^" = (-1)^"^% if 
is odd. In particular, if ~ 0, this factor reflects the 
fact that a hole band is translated by tt in relation to an 
electron band. 

The solution of the above model is a little-trickier than 
that of the usual [/ = oo Hubbard mode£3 due to the 
term crjVe m t^. Its solution is easier to understand if one 
considers first the application of the Hamiltonian in the 
subspace of states with the same configuration of the a- 
spins, \ai, . . . , <jn^)- Then the Hamiltonian can be writ- 
ten in more compact notation, dropping the spin index, 



(36) 



with hopping integrals given as above and Q being the 
cyclic spin permutation operator. 

Consider a general state with no link at site 1, 
|ai, . . . , Cat-; cti, . . . , un^)- If we redefine these states in 
the following way. 



|ai, 



, Cat-; t, • . • , (Tn^) 



(37) 



with oi > 2, the Hamiltonian within the subspace of 
states with the above spin configurations becomes the 
one given by Eq. pq with the following modifications 



ha 

if. 



(38) 



The hoppings across the boundary do a cyclic permu- 
tation of the spin sequence {a} with the above phase 
factor. We wish to construct now the states that remain 
invariant under such a cyclic permutation, that is. 



(39) 



and Tq,^ is the periodicity of the spin configuration and 
ttc labels the different spin configurations. For exam- 
ple, the spin periodicity in | o o | is 3. (f>' will 
be the effective flux felt by the noninteracting fermions. 
This problem is equivalent to solving a one-particle tight- 
binding model for a chain of sites with hopping con- 
stant tj = tg^/'^(^+<^i-<^N^)tP ^ wiih the correspondence 
\i) = Q*-i|{z/}). The total flux through this tight- 
binding chain is 



bi = Tq, — ' —iF. 

^ Ni+ Ni 



(41) 



The solution is obtained after a gauge transformation so 
that tj e*'^^/''°'=t. The gauge transformation depends 
on the ^-spin configuration, but the tight binding eigen- 
values only depend on the total fiux. The eigenstates will 
be Bloch states \ac^ Qc) (m the cyclic permutations) with 
Qc — n{2'K / Va^) , with n = 0, . . . , r^^ — 1. This resolution 
is rather similar to that of the JIubbard model with flux 
which has been treatedJjti_Rcf.llJ. 

Its solution is knowrJlSo and the eigenvalues of Hi for 
L odd are given by 



E{{~k}, (7e, P) - -2t ^ COS (k + a^ + t 

i—i ^ 



with 



a 



Ni+N^ 



(42) 



(43) 



If L is even, there is a tt/L correction in the argument 
of the cosine due to the term (—1)^'*. Note the sign 
change within the cosine argument when compared with 
Eq. This sign change just reflects the "particle-hole" 
transformation which is implicit in the fact that now the 
links (— •) are mapped onto holes. 

Now, the total momentum P has to be determined 
as a function of {fc} and qc- The following condition is 
obtained from the phase acquired by a eigenstate under 
the translation of two real sites or a pseudo-site, 



2P 



p 

a— 
L 



^) (mod 2tt) 



(44) 



which is easy to understand examining the translation 
of a component of the eigenstate which does not have 
pseudo-particles at site L and therefore does not suffer a 
circular permutation of the pseudo-spins. Obviously, the 
components that do not satisfy the previous assumption 
will lead to the same result since the overall eigenstate is 
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invariant by translation. This relation can be written in 
a simpler form 



1j . 1j 

l—l 



(mod 27r) 



(45) 



As in the previous section, the set of pseudo-momenta 
{k} and the pseudo-spin momentum Qc are not enough 
to define totally P. 

The spin-charge factorization of the U — oo Hubbard 
model translates into a decoupling of the degrees of free- 
dom describing the configuration of domains of solitons 
and antisolitons and the degrees of freedom associated to 
the presence of "holes" moving through these domains. 
This factorization and the mapping presented in this pa- 
per are illustrated in Fig. |l]. 



Following the same steps, this leads to the modification 



-<Tii0g-l/2-(l+(Ti-<Tjv~)iA'.^gl/4-(l+<Tjv-)(l+o-i)-i.# 

La 



This phase term generates an extra fiux contribution 
through the {N-^ tight-binding chain which is given 

by 



TT 



a— — h ■ 



(47) 



where N^^i is the number of pairs aa' in the sequence 
of spins obtained with our mapping. For example, in 
Fig. I, N^^ = 1, Nil = 3 and ^ 1. Note that 
Ni = TVtt + Nu, Ni = Nil + ^Ti. ^Ti = ^iT- The 
Hamiltonian is simple to diagonalize and the eigenvalues 
are given by Eq. tt3 with 



A. Flux dependence 

Assume now that the chain is pierced by an external 
flux (/), that is, the Hamiltonian is given by Eq. ^ with 
t te"^l^. This problem can be solved following the 
same procedure as for = with an extra step. This 
step is equivalent to the gauge transformation 

which carries all the phase to hoppings at the boundary, 
tj — > i; j ^ L, th — > te*"^. Let us show how this can be 
done for the mapped Hamiltonian. We modify the state 
invariant by translation in the following way, 



(46) 



Now note the following. 



)p = e^^e-'^= 



(.-..... .)p = e*^e-'^=*e**(. .-.••••) 



• '-jp; 

p- 



Therefore, we will have an extra phase term in the hop- 
pings displayed in the previous tables which involve a 
translation. Furthermore, a hopping of a link (••) at the 
boundary implies a hopping of an electron in the same 
direction while the hopping of a link ( ) implies a hop- 
ping of an electron in the opposite direction. For zero 
external flux, this distinction would be irrelevant, but 
for a finite flux, it leads to a spin dependent phase of 
the hopping integral e~'^i*"^. Following exactly the same 
procedure, we arrive to the same stage of Eq. |3^ with the 
following modiflcations 



'Ll 
til 



g-o-„.-iAr,.£g(l+o-iv,-)-i'>/2^ 



and 



E{{k}) ^ E [ Ik - 13 



(48) 



Ni+Ni 



Ni-Ni 
Ni+Ni 



(49) 



Such expression for the flux dependence should be ex- 
pected since solitons and anti-solitons in the strong cou- 
pling limit can be viewed as hard-core particles with op- 
posite charges and a simple spinless model of hard-core 
particles with opposite charges in a magnetic flux would 
exhibit precisely this flux dependence of the eigenvalues. 
It is easy to show that if iV| = or = 0, /3 = ±1 as it 
should be. One can interpret /3 as the effective charge of 
the carriers. Note that this renormalization of the flux 
dependencepwas also found for the strong coupling Hub- 
bard modelEil with precisely the same form. 



B. Higher order corrections 

The second order corrections can be obtained consid- 
ering virtual hoppings that create or destroy a soliton- 
antisoliton pair. For a given low lying eigenstate with 
Ni — 0, this leads to a energy correction of the form 

— {n^ln^+ll + (1 - n,+ii){cl^c,+2i + c|^2iCa)>- (50) 

When Ni = 0, the energy correction is of the same form. 
If Ni + Ni = L, the second-order corrections can be 
mapped on a Heisenberg spin model. In the general case, 
the energy correction can be written as an average over an 
operator that creates (or destroys) a soliton-antisoliton 
pair and destroys (or creates) also a pair which may or 
may not be the one created (destroyed), leading to long 
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range hopping of these pairs with or without exchange 
of the pair. A closer mapping than that onto the U ~ 
oo Hubbard model is suggested at this level, since the 
above corrections are also present in the charge sector of 
the J7 3> t Hubbard model, if the spin configuration is 
restricted to be Neel like with momentum = 0. In this 
case, long range hopping of a hole- "double occupancy" 
pair is also possible and one may think of doubly occupied 
sites, holes and singly occupied sites (with anqs — Neel 

configuration) as equivalent to (••), ( ) and empty 

sites in our reduced chain. The flux depeadence of the 
eigenvalues also suggests such a picture.EJ We will see 
that such a picture agrees with the transport properties 
of the t-V model. 



IV. COMPARISON WITH BETHE ANSATZ 
RESULTS 

Our results can be li»k|sd,to those obtained with the 
Bethe ansatz technique lla-EJ In the following, we adopt 
the notation of Ref.EJ. The Bethe ansatz solution is 
characterized by a set of bands 7 (with 7 = 0,1,...), 
with non-trivial relations for the total number of avail- 
able "single-particle" states in each band dc.-, and for the 
total number of occupied states in each band, iVf,,^ . The 
energy associated with an occupied state in 7 ^ band is 
of order V and therefore, the 7 = band is the free car- 
rier band obtained in our picture. The high energy bands 
(7 7^ 0) are related to the remaining degrees of freedom 
associated with the possible configurations of the links 

(••) and ( ). This is similar to the strong coupling 

Hubbard model case where the high energy BA bands 
are clearly linked to the_flossible configurations of holes 
and double occupancies. O 

In Table 1 of Ref.ll3, we see that the low lying states 
(A^o 7^ 0, = 0, for 7 > 0) are those of a chain with 
a reduced size L = dc^o — L ~ A^g and number of holes 
given by Nj^ — iV^Q — L — Ne, which agrees with our 



equivalent findings in Section HE . Noting that in Eq. 24 
of Ref.y, 



27r 



27r 



LNe {L - Ne)Ne 



27r 



LiL^NeY 



the eigenvalue expression, Eq. 23 of RefEl, becomes ex- 
actly the same as our Eq. ^ and the same can be said 
for the flux dependence of these eigenvalues. 

The high energy states are more complex since they 
are characterized by a non-zero occupation of the high 
energy bands. Let us assume for simplicity that only one 
of the high energy BA bands (7 7^ 0) is occupied. The 
effective size ofrthe chain and the number of holes in the 
7 = band aretS 



L = 4,0 = L - A^^ + (7 - l)A^c,7, 



(51) 



where in the last equation, we have used the fact the 
holes in the 7 = band are in our picture the particles. 
Relating these two equations to the definition of these 
quantities in our picture, one obtains 



N- 



(52) 



These relations can be confirmed calculating the total 
number of electrons 



A^e = A^cO + (7 + l)^c,7 

= A^_. -N^i+N,,+N^i. 



(53) 



These relations imply that 7 has a very simple physical 
meaning in the strong coupling limit, it is the size of 
the clusters of links (••). Since only one BA band is 
occupied, all clusters have the same size and the total 
number of these clusters is A^| ^ . The total number of links 
(••) is then obviously 7 • Af|t— .|lQiis type of excitations 
form the so-called BA strings^jO and in particular, the 
string associated |-with an occupied state in band 7 has 
length 7 (see Ref.c2 for an explanation of these B A string 
excitations and of the precise meaning of string length). 
We see now that, for t-V model in the strong coupling 
limit, a string is simply a cluster of links (• •) in the 

configuration of links (• •) and ( ). In the general case, 

several BA bands are occupied implying a configuration 

of links (• •) and ( ) where the links (• •) combine into 

clusters of several lengths. 



V. TRANSPORT PROPERTIES 

The transport properties of one-dimensional models 
have acquired a renewed jiiterest recently due to a con- 
jecture by Zotos et ac3~c3 that integrable models with 
zero Drude weight at zero temperature are ideal insula- 
tors, that is, the Drude weight remains zero also at finite 
temperatuiie. Based on qualitative arguments, Zotos and 
PrelovsekE^ have also stated that, in the particular case 
of the strong coupling half-filled t-V model, such temper- 
ature independence is present even for finite size chains. 
This has, been confirmed by a Bethe ansatz study of this 
model.li3 These results can be easily rederived with our 
solution and they are simple consequences of the /3 pref- 
actor in the flux dependence of Eq. For instance, the 
current operator 



J 



can be obtained at finite temperatures from 



Z d{ct)jL) 



(54) 



(55) 



and therefore, if all eigenvalues are flux independent, 
the current will be zero whatever the temperature value. 



8 



Note this is a stronger absence of current than the usual 
situation, which may occur also in metallic systems, 
where the zero average of the current operator results 
from the fact that the positive energy slopes being ex- 
actly compensated tejrtte negative ones. Also, the charge 
stiffness is given byn3E3 

^-^2>^ z di<p./Lr ^^^^ 

and as for the current, if all eigenvalues are flux indepen- 
dent, the Drude weight remains zero at finite tempera- 
tures. A eigenvalue in order to be flux independent must 
have — = 0. It is easy to show that is indeed the 
case for half-filled states. For these states, Ne = Nh and 
since - Ni ^ Ne - Nh, /3 = 0. 

In the following, we illustrate the use of our solution 
with a study of the optical conductivity. The real part 
of the conductivity cr(w) is given by 

Cr{uj) — 2t:Dc6{uj) + areg{'^) (57) 

with 

1 — e^^^ TT v-^ 
o-re3(w) = y > Pn\{n\J\m)\^ 

Lu Li ■^"^^ 

■5{i^ - Era + En) (58) 

where pn is the Boltzmann weight. 

At half-filling, the ground state of the t-V model is in- 
sulating (-Dc(O) — 0) and doubly degenerate, one state 
having momentum 0, the other tt. Both states have 

7V_, = L/2 {L even) and iV,, = N = 0. The current 

operator applied to a ground state induces transitions to 
states with N,, = N = 1 and N_, = L/2-2. 

One should note that when determining the optical 
conductivity at finite temperature, one has to calculate 
matrix elements of the current operator between states 
with iV,, and N^, = N,, -I- 1 in order to obtain the upper 
band part of the optical conductivity. The low frequency 
region is given by matrix elements of the current operator 
between states with the same number of links. Clearly, 
the contribution of the states \n) = |7V„ 7^ 0) will be very 
small as its Boltzmann weight is p„ ~ e~^^"^ /Z and we 
will only consider temperatures T <^ V/ks- So, the sum 
over \n) becomes a sum over all states |A'',, — 0) and the 
sum over |m) becomes a sum over all states |A^,, — 0) for 
the low frequency conductivity and a sum over | A^,, = 1) 
for the upper band part of the conductivity. That is, 
we can write <Treg = Co + fi, where CTq will be the low 
frequency conductivity (w ~ t) 

CTo(t^) = - J Y] PN..=0 (59) 
N',, = a 

N,, = 

\{N:.=0\Jo\N.. = 0)\^S{u;-Em:,=o + En..^o) 
and where ai will be the high frequency conductivity 



f^l(^) = y 2^ PN..=0 (60) 

«;.=! 

N,,=0 

\{N:, = l\Ji\N.. = 0)\H{u; - E^,^^, + En..^o) 

where J — Jq + Ji, Jq being the part of the current 
operator which does not alter the number of links and 
therefore, commutes with the strong coupling Hamilto- 
nian, 

Jo = it'^i^ - rii+2)cJ+iCi(l - "i-i) 

i 

+ it^^ ni+2c\^iCini^i + he, (61) 

i 

and Ji being the sum of terms in the current operator 
which induce transitions between states such that their 
energies differ by V , 

Ji = it'^i^ - nj+2)c|_^iCjni_i 

i 

+ it'^ni+2cl^iCi{l - Tii^i) + he. (62) 

i 

Since Jq commutes with the Hamiltonian, (To(w) = 0. 

In this paper, we will only evaluate ai{uj) at zero tem- 
perature and half-filling {Dc — 0), 

n:.=i 

(1 (<. = 1; P = 0| Ji I A^_. = A^,; P = 0) 1' 

+ = = Ji|A^-. = A^e;P-7r)|') 

• S{uj - £^jv;.=i) 

where |A^_. ^ Ne;P ^ tt) and |A_. ^ N^^P ^ 0) are 
the two possible ground states at half-filling. Well, 

Ji|A_. = AeP = 0> = itVI|l,2;gc(iT) =7r;P = 0), 
Ji|A^_. = AfeP = tt) = ii\/I|l, 2; g,(it) = 0; P = tt), 

so we only need to calculate the overlap of the eigenstates 
with the above states. Note that for these eigenstates, 
Ni ~ N^^ — 1 and therefore N^-^ — 1. Consequently, 
a = 0. The respective eigenvalues are given by 

E{ki,~k2)^V -2tY,cos(k, + '^^^^\ (63) 

[2 L ) \ 2 ) 

These eigenstates are given by 
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From Eq. ^ and noting that 2P = (mod 27r) , one has 
ki + k2 = 2tt/L (mod 27r) ii Qc — and fci + A:2 = (mod 
27r) if Qc = TT. Well, ci has two contributions of the form 



sin2 ^L_^ .5[c.-i?(fci,fc2)], 



(65) 



which taking into account the above conditions can be 
written as sin^(fc) • (5 [w — F + 4i cos(fc)] in both cases. The 
ground state momentum (P = or P = tt) is irrelevant 
as expected. So, 

ai ~ ^ sin^ {k)-S[Lu-V + At cos(fc)] 



1 - 



LU~V 

At 



N 



Id 



LJ-V 



(66) 



where N^'^ is the density of states of a one-dimensional 
tight-binding model. This density of states is non-zero 
between —2t and 2t and has inverse square root diver- 
gences at ±.2t. Therefore, the optical conductivity will 
be characterized by absence of weight between zero and 
V — At, which is the optical gap. At the extremes, the 
optical conductivity goes to zero as ci ^ -\/|cj — y| — At, 
that is. 



V 



At 



-At<Lu-V< At. 



(67) 



This isjjrecisely the dependence obtained by Lyo and 
GalinarLElEj for the stopsc coupling Hubbard model with 
a Neel ground state.c3£d Again, the t-V model seems 
to behave as the strong coupling Hubbard model with a 
fixed Neel spin configuration. The optical conductivity 
of the strong coupling t-V model has also been recently 
studied with arCpnjunction of Bethe ansatz and confor- 
mal invariance,c3 which has allowed the determination of 
the exponent for the frequency dependence immediately 
above the absorption edge. The exponent obtained was 
1/2, in agreement with our results. 



VI. CONCLUSION 

In conclusion, we have presented a non-Bethe-ansatz 
solution for the strong coupling t-V model with twisted 
boundary conditions (or equivalently, the strongly 
anisotropic Heisenberg model). We have found that this 
model can be described in terms of soliton-antisoliton 
configurations and non-interacting particles moving in 
a reduced chain threaded by a fictitious flux generated 
by the previous configuration, but also containing a 
term proportional to the total momentum of the non- 
interacting particles. The flux dependence of the eigen- 
values remains unchanged for the low lying states, but 
is reduced for intermediate energies reflecting the renor- 
malization of the charge of the non-interacting particles. 



Much of the previous picture was obtained with a simple 
mapping of this model onto the U — oo Hubbard model. 
However, the t-V model remains simpler than the U ^ oo 
Hubbard model, since no large spin degeneracy is present 
in the low energy sector. This allows a much easier calcu- 
lation of correlations. As an example, we presented the 
simple calculation of the zero temperature optical con- 
ductivity of this model at half-filling. |_. 

The arguments by Zotos and PrelcHsseka and the 
Bethe ansatz studies of the t-V modellla in the strong 
coupling limit have been confirmed here. All states ob- 
tained from the half-filled insulating ground-state by suc- 
cessive applications of the single particle hopping oper- 
ator have energies independent of the external magnetic 
flux and consequently, the Drude weigth remains zero 
even at finite temperature. However, one should note 
that the results presented here do not exclude a posi- 
tive charge stiffness at finite temperature resulting from 
the t'^/V corrections to the strong coupling Hamiltonian. 
Obviously, the magnitude of the charge stiffness (if non- 
zero) resulting from these corrections will be at the most 
of the order of t'^/V. The flux independence is closely 
linked with the soliton-antisoliton(g) "free particles" fac- 
torization at intermediate energies. It is curious to note 
that such a factorization in the charge sector for interme- 
diate energies, is also present in the strong coupling Hub- 
bard modehliJ The extended Hubbard model can also be 
solved in the strong coupling limit taking a similar path 
to that presented here.Ej 

Finally, note that we have assumed implicitly through- 
out the paper that V was positive, but, obviously, the 
solution is valid for both V/t ^ oo and V/t ^ —oo. For 
t = 0, the eigenstates of the model are the same, inde- 
pendently of the sign of the nearest-neighbor interaction, 
but obviously the respective eigenvalues are symmetric 
for V positive or negative. For t 0, the projection of 
the kinetic energy operator in the degenerate subspaces 
is independent of the sign of the interaction and so will 
be the diagonalization of this operator in each degener- 
ate subspace. Therefore, the eigenstates are the same for 
V/t — > ±oo with the respective eigenvalues being given 
by Eq. 42 plus or minus N„ • |T^| , according to the sign of 
the interaction V. In particular, at half-filling, the phase 
separated state (which is the highest energy state when 
V is positive) and the two charge ordered groundstates 
trade places when V is negative, i.e., the phase separated 
state becomes the ground state and the two charge or- 
dered states become the highest energy states. 
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FIG. 1. The mapping is illustrated in this Figure. In the 
first chain, the circles stand for occupied sites while the small 
dots stand for empty ones. The introduction of the near- 
est-neighbor interaction leads to a further factorization of the 
wave function describing the charge degrees of freedom. 
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